Anomalous Scaling in a Model of Hydrodynamic Turbulence with a Small Parameter 



Daniela Pierotti, Victor S. L'vov, Anna Pomyalov and Itamar Procaccia 
Department of Chemical Physics, The Weizmann Institute of Science 
Rehovot, 76100, Israel 



OS 
On 
On 



(N 



> 

o : 

m . 

o ■ 

m ■ 

o : 

Os . 
Os ■ 

: 

>v 

■ 

o : 

O ■ 



The major difficulty in developing theories for anomalous scaling in hydrodynamic turbulence is 
the lack of a small parameter. In this Letter we introduce a shell model of turbulence that exhibits 
anomalous scaling with a tunable small parameter. The small parameter e represents the ratio 
between deterministic and random components in the coupling between N identical copies of the 
turbulent field. We show that in the limit N — » oo anomalous scaling sets in proportional to e 4 . 
Moreover we give strong evidences that the birth of anomalous scaling appears at a finite critical e, 
being e c pb 0.6. 



The statistics of the small scale structure of turbulence 
is characterized by "anomalous scaling" meaning that 
correlation functions and structure functions of velocity 
differences across a scale R exhibit a power law behavior 
with scaling exponents that are not correctly predicted by 
dimensional analysis. In the last few decades there have 
been many attempts to compute the scaling exponents 
of turbulent fields from the equations of motion. In the 
context of simplified models of passive scalar advection 
it was discovered that there exist natural small parame- 
ters that allow direct computations of anomalous scaling 
exponents 0-^]. In Navier Stokes turbulence and also in 
simplified models like shell models the progress in com- 
puting scaling exponents was slowed down by the lack of 
a small parameter. It is thus worthwhile to consider mod- 
els of turbulent velocity fields in which a tunable small 
parameter can be introduced and used to advantage. 

We propose to introduce a small parameter via the 
coupling between N copies of a field u n which satisfies 
the dynamics of the Sabra shell model of turbulence in- 
troduced in 01: 



du n (t) 
dt 



ak n+ i u* +1 u„ +2 + bk n u*_ x u n+ i 
ck n — i 

U n -2«7i-X - U n + f n (t) . (1) 



Shell models are simplified dynamical systems con- 
structed such that the complex number u n represents the 
amplitude associated with the Fourier transform of the 
velocity field u(r) with "wave- vector" k n . Rather than 
considering the full k space and all the nonlinear inter- 
actions one allows for only one-dimensional k vectors on 
shells spaced such that k n = koX n , with A being the spac- 
ing parameter, and local interactions. In Eq. (^) v is the 
"viscosity" and /„ (t) a random Gaussian force restricted 
to the lowest shells. The parameters a, b and c are re- 
stricted by the requirement a+b+c — which guarantees 
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the conservation of the "energy" 
the in- viscid, unforced limit. 

We now want to generalize this model to one which 
consist of N suitably coupled copies of it. In order to 
do that we need to consider separately the equations for 
the real and imaginary parts of u n . This procedure guar- 



antees that the obtained model converges to the original 
Sabra model in the limit N — > 1. The copies are indexed 
by i, j or £, and these indices take on values — J, . . . , + J, 
2 J + 1 = N. The ith copy of the velocity field is de- 
noted .a. In this notation a = ±1 refers to the 

real and imaginary parts of u n respectively. Let D™*' be 
the coupling between copies, which will be chosen later. 
Equations (Fy) for a collection of copies are 



du 



it 



dt 
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c > n 1 n— 2,ct' n— 1,(7 



where 
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To proceed we note that the index t is defined mod- 
ulo N, and introduce a Fourier transform in the "copy" 
space, defining the collective variables: 
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N ^u^e X p { N 



(5) 

l=-j 

Note that the index a is also defined modulo N = 2J+1. 
It is convenient to consider values a within "the first Bril- 
louin zone" , i.e from — J to J. We will refer to it as the 
a-momentum. Since tin] CT is real, u^" = u"^. In "cv- 
Fourier space" Eqs. (||) read 

\-AJ yj-l .y-j ^ 

/3,7 



dt 



V < <,* + fn,a ■ (6) 



n-2,a' w n- 



1 



where A a .p is the Kronecker symbol. Observe that we 
use Greek indices for components in a-Fourier space, and 
Latin indices for copies in the copy space. As a conse- 
quence of the discrete translation symmetry of the copy 
index [i] Eqs. (|^) conserve a-momentum modulo N at the 
nonlinear vertex, as one can see explicitly in the above 
equation. The coupling amplitudes $ a '' 9 'T in these equa- 
tions are the Fourier transforms of the amplitudes D^^. 
We choose the coupling amplitudes according to 



a,/3,7 



= -$=[e+y/l-e* ^> a ^} , 



(7) 



where ^ a ^-^ are quenched random phases, uniformly 
and independently distributed with zero average (cf. 
the "Random Coupling Model" (RCM) for the Navier- 
Stokes statistics J5| and the identical symmetry condi- 
tions there). Consequently for e = the model reduces to 
the RCM and exhibits normal scaling (K41) for N — > oo. 
It was in fact understood |||(| that in the limit N — > oo 
the direct interaction approximation (DIA) becomes the 
exact solution of the RCM. Moreover a proper analysis of 
the DIA approximation leads to normal scaling for those 
systems in which sweeping effects are removed or absent 
by construction like in shell models 

On the other hand for e = 1 the coupling coefficients 
in the a- Fourier space (Q) are index- independent. This 
corresponds to uncoupled Eqs. (||) in the copy space, 
because in this case — SijSij. Thus for e = 1 we 

recover the original Sabra model with anomalous scaling 
j|. Our choice of couplings (]?]) allows an interpolation 
between the normal K41 scaling for e = (at N — > oo) 
and the full anomalous scaling of Sabra model for e = 1 . 
A model of this type was proposed in the context of 
Navier-Stokes statistics by Kraichnan in || and analyzed 
by Eyink |6| in terms of perturbative expansions. 

Our aim in this Letter is to present numerical re- 
sults which show that for small values of e the model 
exhibits small anomalous corrections to normal scaling 
and to present theoretical arguments to rationalize the 
functional dependence of the anomalous corrections on 
e. We measured the scaling exponents of the structure 
functions: 



S2 P (k n ) 



i— — J a 



(8) 



The exponents C%p have been calculated by a linear fit in 
the two decades inertial range, see Fig. [l| The equations 
of motion (||) with 28 shells, a = 1, b = c = —0.5, were 
integrated with the slaved Adams-Bashforth algorithm, 
viscosity v = 4 x 10~ 9 , a time-step At = 10 -5 . The forc- 
ing was subjected on the first two shells, chosen random 
Gaussian with zero average and with variances such that 
02/ci = 0.7 (in order to minimize the input of helicity 
Pi which leads to period two oscillation in the structure 
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FIG. 1. Log-log plot of the structure functions S p (k n ) vs 
k n for p=2,3,4, e = 0.8 and N = 25. 
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FIG. 2. 5^2 = <2 - 2/3 vs 1/N for e = (circles), e = 0.8 
(squares) and e = 1 (diamonds) for N from 5 to 25. The point 
at 1/N = for the e = curve is the theoretical prediction of 
the RCM for N -> 00. 

functions). Averages were taken for a time equal to 250 
eddy turnover times for the case N = 1. The averag- 
ing times were decreased when the number of copies in- 
creased, taking into account the faster convergence times 
in these cases. The quality of the scaling behavior and 
of the fits is demonstrated in Fig. [j]. 

To substantiate the birth of anomalous scaling at e > 
we simulated the model for different values of e and of N. 
We are interested in the values of the scaling exponents 
for very large values of N (00 in theory). In Fig. || one 
can see the plot of the value of the anomalous corrections 
to Kolmogorov scaling, 8(2 = (2 — 2/3, as function of 1/N 
for e = 0.8 together with the same curve for e = and 
for e = 1 for N ranging from 5 to 25. While for e = the 
corrections to Kolmogorov scaling go to zero, for e = 0.8 
and for e = 1 the corrections converge to a finite value 
which increases with e. 
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FIG. 3. 6Q2 — C2 — 2/3 vs e 4 (circles), vs e 2 (squares) and 
vs e (diamonds) together with linear, quadratic and quartic 
fits respectively for N=25. 
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FIG. 4. S( 2 = C2 — 2/3 (circles), 8(4 (squares) and 8(e (di- 
amonds) vs e 4 for N — 25. 

The random phases in the couplings were chosen with 
respect to a uniform probability with zero-mean at the 
beginning of each simulation. The rigorous procedure 
for quenched disorder would call for taking averages over 
different runs with different choices of the couplings. We 
did not do that, but rather checked that self-averaging 
is already valid for N = 5, for e = 0.8 (small random 
component in the couplings) and within our numerical 
precision. For e = self-averaging occurs only for large 
numbers of copies. 

The most interesting aspects of the numerical findings 
are the dependence of the anomalies on e and the question 
whether anomaly appears for any e > or only above a 
critical value e c . The first issue is settled with sufficient 
clarity in Fig. in which we show the behavior of 
(for N=25) as a function of e, e 2 and e 4 with the re- 
spective linear, quadratic and quartic fits. One can see 
clearly that the anomalous corrections SC2 go to zero like 
e 4 . The same behavior is exhibited by 5(^4 and 8<^e as one 
can see in Fig EL 
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FIG. 5. 5^6 = C6 — 2 vs e for N = 25 with the maximal 
and minimal slope lines obtained for the linear fits in the two 
ranges e 4 G [0.13, 0.88] and e 4 g [0.23, 0.88] . All the values of 
the anomalous corrections have been shifted by the average 
of the first five points. 

The existence of a critical value e c is harder to set- 
tle. In the following we show arguments in favor of the 
existence of a finite e c . First one should notice in the 
magnification in Fig. || that at e = we do not get the 
K41 values S( n = as theoretically predicted. This is a 
result of the finiteness of the number of copies (N = 25 
in Fig. |). 

To establish the existence of a finite e c we proceeded 
as follows: i) To take this into account the finite size ef- 
fects we subtracted the e — value of 5( n from all the 
values of £„ . ii) We calculated the maximal and minimal 
slope line, that is the best fit slope plus and minus the 
standard deviation, for 5C, n vs e 4 in two different ranges: 
e 4 G [0.13,0.88] and e 4 G [0.23,0.88]: that is, we fit- 
ted with the 2 parameter function /(e 4 ) = a n + b n e 4 . 
iii) We fitted the plot of S( n vs e with the 4 parameter 
function 17(e) = a n e + b n e 2 + c„e 3 + d n e 4 in the range 
e G [0, 0.97]. Note that with the "calibration" of the zero 
that we performed 5( n — for e = 0. The \ 2 test is 
much better for the first procedure than for the second, 
where the minimal % 2 is tripled, although in principle it 
should be easier to fit a function with a larger number of 
parameters. 



5Cn 


range of e 4 


slope xlO 3 


e c 


5(2 


[0.13,0.88] 


48± 3 


0.60 ± 0.07 


<5C2 


[0.23,0.88] 


50± 3 


0.62 ± 0.06 




[0.13,0.88] 


-74 ±8 


0.69 ± 0.07 


5(4 


[0.23,0.88] 


-82 ±9 


0.71 ± 0.06 




[0.13,0.88] 


-341 ± 8 


0.59 ± 0.04 




[0.23,0.88] 


-335 ± 6 


0.57 ± 0.05 



TABLE I. Results of the linear fits on different ranges 
for different scaling exponents. The values of e c have been 
calculated as explained in the body of the text. 
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We interpret this result as an evidence for the existence 
of a finite e c . 

In order to have a better estimate of the "zero" level of 
S(n's, and so a better estimate of e c , instead of subtract- 
ing the value of the anomalous correction at e = we 
subtracted the average value of S( n calculated by using 
the first 5 points belonging to the flat region. The value 
of e c with its error has been found by looking at the in- 
tersections of the minimal and maximal slope lines with 
the zero line: the line of the value of S( n for e = with 
its error (see Fig. ||). Table | exhibits the resulting values 
for the various exponents and their linear fits. Note that 
all the values of e c coincide within the error bars. 

To understand the linear dependence on e 4 we turn to 
the exact equations that are satisfied by the n-order cor- 
relation functions, which are defined as suitable averages 
over time and over all the replica. The second and third 
order correlation functions are defined as: 

ftfe;t-t')4D^W)) ) (9) 

a,cr 

Fs (k n ;t,t',t") = ± J2 E $Q "' Q ' Q ( 10 ) 

a, a' ,a" cr,<r' ,cr" 

«-!,.' (*'K+V(*")) • 

and similarly for higher order correlation functions and 
for the response functions (for details see ||). 

The equations of motion for these objects can be writ- 
ten down explicitly. For example, 

d 

-^F 2 {k n ,t) = j atn+1 F 3 (k n +i;0,t,t) +fb, n F 3 (k n - 1 t,0,t) 

+7c,n-iF 3 {k n -i; t,t,0), (11) 
d 

-^-—F 3 {k n ,ti,t2, t 3 ) = r y a ,nF4(k n , fcr»+l> kn, fcri+li *1) *1) *2> *3) 

at i 

+lb,n-iFi(k n ^2,k ni k n , k n+ i, t\, t\, t 2 , t s ) 

"HTc.n— 2F4 {k n — 3, k n _2 > km i 

In order to close such equations and to attempt to solve 
them, one needs to express a higher order correlation 
function in terms of lower order statistical objects. In ||] 
it was shown that this can be done in the present context 
in a controlled fashion. In other words, it is possible to 
express, say, F4 in terms of F2, F$ and second and third 
order response functions. In doing so one of course leaves 



out information about F4 that cannot be possibly repre- 
sented in terms of lower order objects. Yet, the main 
result of the analysis of || is that the neglected terms in 
this procedure are of 0(e 6 ) whereas the retained terms 
are of 0(1) and of 0(e 4 )\ Since we understand that for 
e = the anomaly must vanish, we expect the anomalies 
to be proportional to e 4 . We interpret therefore the nu- 
merical results shown in Fig || and Fig ^ as an excellent 
confirmation of this theoretical expectation. 

To summarize: we introduced a shell model of hydro- 
dynamic turbulence in which a small tunable parameter 
e exists. We know from theoretical and numerical results 
that the model exhibits normal scaling for e = 0. The 
most important results of this letter are the existence of 
a finite e c at which anomalous scaling appears and the 
behavior as e 4 of the anomalous corrections for e > e c . 
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